Wave localization in binary isotopically disordered one-dimensional harmonic chains 
with impurities having arbitrary cross section and concentration 



^1" 
O 

o 

(N 



CZ5 



K.A. Snyder 

Materials and Construction Division, National Institute of Standards and Technology 

T.R. Kirkpatrick 

Institute for Physical Science and Technology and Department of Physics, University of Maryland 

(Dated: February 2, 2008) 

The localization length for isotopically disordered harmonic one- dimensional chains is calculated 
for arbitrary impurity concentration and scattering cross section. The localization length depends 
on the scattering cross section of a single scatterer, which is calculated for a discrete chain having 
a wavelength dependent pulse propagation speed. For binary isotopically disordered systems com- 
posed of many scatterers, the localization length decreases with increasing impurity concentration, 
reaching a mimimum before diverging toward infinity as the impurity concentration approaches a 
value of one. The concentration dependence of the localization length over the entire impurity 
concentration range is approximated accurately by the sum of the behavior at each limiting concen- 
tration. Simultaneous measurements of Lyapunov exponent statistics indicate practical limits for 
the minimum system length and the number of scatterers to achieve representative ensemble aver- 
ages. Results are discussed in the context of future investigations of the time-dependent behavior 
of disordered anharmonic chains. 

PACS numbers: 45.05.+X, 45.30.+S, 46.40.Cd, 62.40.+d 
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I. INTRODUCTION 

The length scale over which Anderson localization^ oc- 
curs in harmonic disordered chains can be expressed as a 
function of the ensemble-averaged system resistivity^. 
Using scaling arguments^ the ensemble-averaged resis- 
tivity of harmonic systems having binary isotopic disor- 
der (single- valued impurities) can be expressed as a func- 
tion of the scattering cross section of a single impurity 
and the impurity concentration. Binary isotopically dis- 
odered harmonic chains also have the property that in the 
limit that the impurity concentration approaches unity, 
the system is once again "pure," and the Anderson local- 
ization length diverges toward infinity. As a result, the 
Anderson localization length passes through a minimum 
at intermediate concentrations-. These properties of a 
binary isotopically disordered harmonic chain (BIDHC) 
make these systems both tractable and interesting. 

In this work, an approximation is developed for the 
localization length of a BIDHC with arbitrary impurity 
cross section and concentration. In addition, the result 
incorporates the discrete nature of a classical chain. Sim- 
ilar studies have been performed previously, but with im- 
portant differences. Bourbonnais and Maynard^ studied 
energy transport in one-dimensional systems having iso- 
topic disorder, but the impurity masses were not single- 
valued. The results of Azbel and Soven& on quantum 
systems having binary isotopic disorder were based on 
a continuum solution for impurity locations constrained 
to exist on lattice sites. Although the Azbel and Soven 
result applied to short wavelengths that may exist be- 
tween the impurities, the results do not incorporate the 
additional features of a discrete mechanical system. 

This study serves as an introduction to future work 



on energy transport in binary isotopically disordered 
anharmonic chains. These anharmonic chains serve as 
a useful model for, among other things, line width 
broadening 7,8,9 . To perform that work, numerical inte- 
gration will used to study the time-dependent nature of 
these systems. To be practical, the initial conditions will 
require sufficiently short wavelengths and high impurity 
concentrations to keep integration times manageable. In 
addition, the results of the study will incorporate the dis- 
crete nature of the mechanic chains so that these effects 
can accounted for in the results. 

In this study, the energy localization in a BIDHC is 
studied for arbitrary displacement wavelength, impurity 
concentration and scattering cross section. Disorder is 
effected by changing randomly selected masses by a fixed 
amount. A continuum Kronig-Pcnney (KP) modeUfi is 
used to develop a general expression for scattering cross 
section, and the continuum impurity impedance is cor- 
rected for wavelength dependent pulse propagation speed 
in discrete systems. The resulting expression is veri- 
fied by direct numerical integration. The localization 
length of systems with strong scatterers is calculated 
using both the continuum KP model and the MacKin- 
non and Kramer 11 (MK) method. The distribution of 
Lyapunov exponents is studied using the continuum KP 
model, and the minimum requirements are found for sys- 
tem length and number of scatterers to achieve proper 
scaling statistics. The localization length concentration 
dependence is studied using the MK method, and an ana- 
lytical expression is found for its behavior. For arbitrary 
impurity concentration and impurity cross section, the 
localization length in a BIDHC can be approximated by 
invoking a simple ansatz based on an analogy to electri- 
cal systems. The result is accurate for systems having 
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displacement wavelengths at least four lattice spacings 
long. 



II. MODEL SYSTEM 

The physical model used here is the harmonic version 
of the Fermi-Pasta-Ulam (FPU) 1 — chain that is composed 
of discrete springs and masses. The masses mi are on a 
lattice with spacing a and interact via nearest neighbor 
springs with force constant K. Disorder is effected by 
changing the background mass to by a fixed amount m + 
with probability c. To simplify the results, all lengths are 
scaled by the lattice spacing a. 

For a system composed of N masses, each character- 
ized by a displacement Xi about the equilibrium location 
and a momentum pj, the Hamiltonian is separable: 
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The real space equation of motion is 

m. 
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The Fourier transform leads to the corresponding equa- 
tion of motion for the energy eigenstate amplitudes Ui(to): 



Ui+l = 
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This is the corresponding Anderson tight-binding model 
for the chain. 



A. Discrete Analysis 

The time-dependent properties of the system were de- 
termined by numerical time integration of Eqn. fusing a 
fourth-order symplectic integrator algorithm (SIA4) for 
separable Hamiltonians. The coefficients were taken from 
Candy and Rozmusi^, and the time step was 1 /200 of the 
natural period. 

The localization length was calculated from Eqn. fus- 
ing the method of MacKinnon and Kramer (MK)Ai. This 
method exploits the statistical properties of the u, so that 
periodic rescaling can be used to improve overall statis- 
tics. 



B. Continuum Analysis 

A Kronig-Penney modeliS is used to develop an expres- 
sion for the scattering cross section of an impurity and to 
study the statistics of the scaling parameter 3 . The con- 
tinuum system analogous to the discrete chain is a ho- 
mogeneous elastic rod having mass density and Youngs 



modulus E. In the absence of impurities, a longitudinal 
displacement amplitude ip(x,t\u>) with angular velocity 
uj will propa gate down the rod with longitudinal velocity 

Q = y/E/ll. 

A harmonic oscillator impurity located at x' will give 
rise to a reactive force due to the impurity impedance Z 
for a wave with angular velocity u: 



dt 2 



E 



dx 2 



-Z(u>) 8{x — x ) 



i/>(x,t\u) (4) 



We assume the solution has a time dependence given by 
an exponential of angular velocity u> (ip = 4>(x\uj) e~ luJt ) 
to obtain 



dx 2 



— 1U! 

— Z{u)5{x- 



4>(x\w) (5) 



where n = to/ci. Equation[S]can be solved analytically for 
a single scatterer or can be solved numerically for a par- 
ticular system composed of a number of scatterers. From 
this numerical solution, one can determine the system re- 
sistivity, which can be used to calculate the localization 
length. 



III. CROSS SECTION 

The scattering cross section a for a single impurity can 
be calculated from the continuum system of Eqn. The 
Green's function for the 1-D Helmholtz equation^ can 
be used to solve for the admittance <fi(x): 
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(6) 



The scattering cross section in 1-D is equivalent to the 
reflection probability: 
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(7) 



To apply this equation to the discrete chain, it must 
be converted from a continuum description to the cor- 
responding discrete description. 

For an FPU system having masses spaced a distance a 
apart, the continuum coefficients can be expressed as a 
function of the discrete properties in the limit a — > Oi^: 



a 



E = Ka 



(8) 



The scattering cross section can now be expressed as a 
function of the discrete system components: 
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\Z(lu)\ 2 + AKm 



(9) 



The impedance Z(u>) of a mass impurity along a one- 
dimensional chain is proportional to the mass m + that is 
added to the background mass mj&: 



Z(io) — — iwTO + /c s 



(10) 
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FIG. 1: Schematic of cross section numerical experiment; each 
line represents the state of the system at the time of mea- 
surement. The time t, shown along left side, is expressed as a 
function of the system length L and pulse propagation speed 
c s ■ The large dots denote the location of two impurities. 



This equation has been modified from the continuum re- 
lation to account for the properties of a discrete one- 
dimensional chain. A displacement wavelength has a 
corresponding wave number k = 2ir/\. For the discrete 
chain, the angular velocity to and relative pulse propa- 
gation speed c s can be expressed as a function of the 
displacement wave number Atii: 



2sin(>/2) 



c s = cos(fe/2) 



(11) 



Finally, substituting Eqn.^]into Eqn. Ogives the scatter- 
ing cross section of a mass impurity in a one-dimensional 
chain: 



(m+ lu/c s ) 2 



(m + uj/c s ) 2 + AKm 



(12) 



A numerical experiment was performed to test the ap- 
plicability of Eqn. E| to FPU systems, and a schematic 
of the experiment is shown in Fig. ^ The system had 
fixed ends and length L. One impurity was located at 
I//3, and another impurity was located at 2L/3. An ini- 
tial displacement was made with wavelength A and total 
length 16A, and the hyperbolic tangent function was used 
to taper the oscillation amplitude from A to zero. Ini- 
tially, the pulse was located at one end of the system. It 
had length 8A, amplitude 2A, and zero initial velocity. 
The time-dependent behavior was determined by numer- 
ical integration using the SIA4 algorithm. 

The measured scattering cross section was determined 
from the energy located in the three regions separated by 
the impurities. The pulse energy located in each three 
regions, denoted by E\, E2, and £3, was calculated at 
four separate times: t = 0, L/6c s , L/2c s , and bL/Qc s . 
At these times, the pulses are located near the center of 
a region. Although redundant with respect to the initial 
total energy, determining the energy at t — L/6c s pro- 
vided a consistency check. In each case, the difference 
between E%(t = 0) and E\(t = L/6c s ) was less than one 
part in 10 4 . 
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FIG. 2: Cross section a as a function of m+Lu/c s for different 
wavelengths A. The solid curve is the analytical result in 
Eqn. 1121 The inset shows additional data near zero. 



The scattering cross section in FPU systems is the re- 
flection coefficient. In this experiment, the reflection co- 
efficient R is calculated from the ratio of energies in the 
first two intervals after one scattering event: 



R 



£i(£/2c s 
Ex(0) 



(13) 



Although the transmission coefficient T could have been 
determined from E2(L/2c s ), a second scatterer was used 
as a more rigorous test of the experiment design and nu- 
merical integrator. The transmission coefficient T for a 
single scatterer was calculated using the energy after two 
scattering events: 



T 



t E 3 {5L/6c s ) 
£i(0) 



(14) 



In all cases, the magnitude of 1 — R — T was less than 
lfj- 3 . 

A comparison of the estimated cross section a in 
Eqn. E| to the measured reflection coefficient R in 
Eqn. El is shown in Fig. [21 for different displacement 
wavelengths A. The results demonstrate that Eqn. IT51 
is an accurate estimate for the scattering cross section a 
for displacement wavelengths as short as 4. Moreover, 
the symmetry about zero for scattering cross section for 
negative values of m + is shown in the inset of Fig. El 



IV. LOCALIZATION LENGTH 

The resistivity pn of a system composed of N scatter- 
ers is the ratio of the system reflection coefficient Rjy to 
the system transmission coefficient T]y2£: 



Pn 



Rn 



(15) 



The localization length of the system is defined in an 
averaged sense. To achieve proper scaling, however, one 
defines its inverse, the Lyapunov exponent 7: 



7 



ln(l + pjv) 



(16) 



More specifically, there is a distribution of values from an 
ensemble of systems. If the systems composing the en- 
semble are sufficiently large, the distribution of 7 values 
will be normal. The localization length £ for a system 
having length L and N scatterers is defined from the en- 
semble averaged Lyapunov exponent: 



r 1 = (7) 



(In (1 + PAf )) 
L 



(17) 



Unless otherwise noted, the symbol 7 shall imply the 
ensemble averaged quantity. 



A. Strong Scatterers 

To perform numerical experiments on anharmonic sys- 
tems of manageable length, the scatters will need to be 
relatively strong. Therefore, weak scattering results will 
not be applicable. Moreover, a means is needed to pre- 
dict the averaged localization behavior of a system using 
only single scatterer information. Because the impurities 
are identical, the scaling law£ can be exploited to express 
system behavior as a function of the resistivity of a single 
scatterer p: 



(ln(l + ptf)) =JVln(l+p) 



(18) 



Substituting from Eqn. El above yields an unbiased es- 
timate for the ensemble averaged localization length (di- 
lute limite) as a function of impurity concentration c: 



(19) 



0), one recovers 
where A is the 



In the limit of weak scattering (c, p 
the expected result = ca = A 
classical mean free path. 

To demonstrate both the effect of strong scatterers 
(cr — ► 1) and the accuracy of Eqn. the ensemble av- 
erage (ln(l + pn)) was calculated from 10 000 systems, 
each having length L and N impurities. The displace- 
ment wavelength was 32, and the scattering cross section 
of each impurity was a. The results from the calculation 
are shown in Fig.[3]as a function of the impurity cross sec- 
tion a. In the figure, the error bars represent the mean 
standard deviation, and many of the symbols lie upon 
one another. As can be seen, all the systems have the 
ln(l + p) dependence that deviates from the mean free 
path approximation {S,' 1 = ca) for cross sections greater 
than approximately 0.2 for all combinations of system 
size and number of scatterers. 
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FIG. 3: Localization length £ , scaled by impurity concen- 
tration c, as a function of impurity scattering cross section 
a for systems having length L and N impurities. The solid 
line is ln(l + p) and the dashed line represents the mean free 
path estimate. The error bars represent the mean standard 
deviation. (Many of the symbols lie upon one another and 
the error bars are typically smaller than the symbols.) 



B. Statistics 

The results shown in Fig. [3] demonstrate that, in the 
mean, systems having a finite density of scatterers have 
the expected behavior. Recent results suggest that sys- 
tems having relatively few scatterers do not exhibit Gaus- 
sian behavior and, therefore, not obey scaling laws. For 
systems having sufficient length and number of scatter- 
ers, the population of Lyapunov exponents is normally 
distributed^, with variance s£ 



2 

1? 



(20) 



Returning to the data of Fig.0 the population of Lya- 
punov exponents 7 was compared to the expectations of 
Eqn. [3U| To assess the "normality" of the data, the in- 
tervals, both above and below the mean, having coverage 
factors^ corresponding to one and two standard devia- 
tions were determined from the population. In addition, 
the population standard deviation A so that it could be 
compared to both its estimated value in Eqn. |20 an d to 
the corresponding coverage interval. 

The results of the Lyapunov exponent statistics calcu- 
lation, from an ensemble with a population of 10 000, are 
shown in Fig. 0] for systems having length L = 16 384, 
N = 128 impurities, and displacement wavelength A = 
32. (Results for other systems were similar, and are omit- 
ted for brevity.) The data are shown as a function of the 
single impurity cross section a. The filled symbols are the 
average value, and the solid line is the estimated average 
value. The error bars represent the intervals that have 
the same coverage factor as one and two standard devi- 
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FIG. 4: Lyapunov exponent 7 statistics as a function of scat- 
tering cross section a for a system with L — 16 384, N = 128, 
and A = 32. Error bars represent coverage factors correspond- 
ing to one and two standard deviations of a normal distribu- 
tion. The predicted intervals ±ls 7 and ±2s 7 (Eqn. 12771 are 
shown, along with observed population standard deviation A. 



FIG. 5: Lyapunov exponent coefficient of variation A/7 as 
a function of the ratio of system length L to the localization 
length £ for systems having different number of scatterers N. 
The solid curve is Eqn. I2l| 



V. CONCENTRATED IMPURITIES 



ations in a normal distribution. The two pairs of dashed 
lines are the estimated standard deviations from Eon. 1201 
The dotted line is the population standard deviation A. 

For small scattering cross sections, the distribution of 
7 is asymmetric, with zero as a lower bound for the cov- 
erage intervals. Interestingly, only the outer intervals 
are asymmetric about the mean. The inner intervals are 
nearly symmetric about the mean, and they have a value 
nearly equal to the population standard deviation. 

As the scattering cross section increases, the results 
begin to exhibit Gaussian behavior. Above a 0.2 scatter- 
ing cross section, the measured intervals, the predicted 
intervals s~, and the population standard deviation A 
all agree. Therefore, agreement between the population 
standard deviation A and the estimated standard devi- 
ation s 7 is as much as measure of "normality" as is a 
careful analysis of the population coverage intervals. 

Exploiting this relationship, the population standard 
deviation A for all the data from the Lyapunov exponent 
statistics experiment are shown in Fig. [S] as a function of 
the system length. For normally distributed populations, 
the estimated coefficient of variation can be determined 
from Eqn. 1171 and Eqn. [23 and has the following conve- 
nient form: 



(7) 



(21) 



This equation is shown as a solid line in Fig. One 
can conclude from the figure that in order for the ob- 
served Lyapunov exponents to be normally distributed, 
the system length must be at least 10 times the localiza- 
tion length, and the number of scatterers must be greater 
than approximately 32. 



For dilute impurity concentrations, the localization 
length decreases with increasing impurity concentration. 
As the impurity concentration c approaches unity, how- 
ever, the FPU system will become a pure system com- 
posed entirely of masses m + to + . At c — 1, the sys- 
tem is once again devoid of impurity and the localization 
length goes to infinity. At intermediate impurity concen- 
trations, the localization length passes through a mini- 
mum. Therefore, the behavior of a system at arbitrary 
impurity concentration cannot be fully characterized by 
the relation in Eqn. EH 

For dilute systems, the localization length £ c ^o is as 
before: 



Cio=C 1 =c<ln(l + p)) 



(22) 



At high concentrations, the localization length £ c ^i has a 
concentration dependence that is proportional to (1 — c): 



(l- C )(ln(l+p')> 



(23) 



The adjusted resistivity p' is for a system having back- 
ground mass TO + m + and impurities with mass m . By 
the nature of the solution using the MK method, and 
given that the systems are harmonic, the frequency u> is 
the same for both systems, but the wavelength and the 
corresponding pulse speed are different: 



k = k \J TO + TO_|_ 

c' s = cos(fc'/2) 

j ^2 



(-m+uj/c' s y 
4K(m + to + ) 



(24) 
(25) 

(26) 



The behavior of the system for all values of impurity 
concentration is conjectured from the electrical analogy: 
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FIG. 6: Localization length £ as a function of impurity con- 
centration c in a discrete system having scatterers with cross 
section 0.2. Solid lines is Eqn. 1271 dashed line is locus of min- 
ima, and dotted line denotes equal contribution from £ c ^o 
and fc-fi. Inset shows same data as a function of (1-c). 



as Lyapunov exponent is to resistivity, localization length 
is to conductivity. If one assumes that at some intermedi- 
ate concentration the behavior is simultaneously express- 
ing itself as two systems with localization lengths £ c _^o 
and £ c _ >!, these two systems should contribute indepen- 
dently to the overall behavior. By analogy to conductors, 
the total localization length £ would be additive: 



(27) 



This equation represents a more complete estimate for 
the localization length that is valid for both strong and 
concentrated scatterers. 



A. Positive m+ 

The accuracy of this approximation is shown in Fig. 
for systems having impurity scattering cross section a = 
0.20 and displacement wavelengths ranging from 8 to 128. 
The localization length £ for these systems was calculated 
using the MK solution. As can be seen, the approxima- 
tion in Eqn. [Z7\is reasonably accurate. The inset shows 
the same data, plotted as a function of (1 — c), highlight- 
ing the separate behavior near c — * 1 for the different 
wavelengths. 

Also shown in Fig. [5] are a dashed line and a dotted 
line. The dashed line is the locus of localization length 
minima as a function of displacement wavelength. Not 
only does the minimum localization length increase with 
increasing wavelength, the concentration at which this 
happens decreases with increasing wavelength. The dot- 
ted line in Figgis the locus of points where £ c _*o = £c->i- 
This locus of points is meant to delineate the cross-over 
point as the system passes from one dominate phase to 
the other. The cross-over point has a stronger concen- 
tration dependence than the locus of minima. 



FIG. 7: Localization length £ as a function of impurity con- 
centration c for m+ = ± 0.8. Inset shows same data as a 
function of (1-c). 



B. Negative m+ 

At a given oscillation wavelength and frequency, the 
impurity scattering cross section depends upon the mag- 
nitude of the change in mass m+. For small impurity 
concentrations, the localization length for a system with 
impurity mass m + ro+ will be equal to one with impu- 
rity mass m — m+. At higher concentrations, however, 
the behavior of systems will differ. 

As an example, the localization concentration depen- 
dence was calculated for two systems with displacement 
wavelength A = 32. The added impurity masses were 
+0.8 and —0.8, and the localization length was deter- 
mined by the MK method. The results of the calculation 
are shown in Fig. das a function of the impurity concen- 
tration c. The estimate from Eqn. I27| is shown as a solid 
line. Within the inset are the data plotted as a function 
of (1 — c). As expected, the behavior of the two systems 
diverge at higher concentrations. 



C. Azbel and Soven Comparison 

To better judge performance, Eqn. 1271 is compared to 
the more rigorous result of Azbel and SovenS, (AS). The 
AS model contains quantum particles interacting with 
delta function potentials that have strength V and are 
located at random integer values of x. For the AS sys- 
tems, the value of V is equal to -1, and the scattering 
cross section <tas of each scatterer is a function of the 
particle wave number k: 



1 



CAS = 



1 + 16sin 2 (fc/2) 



(28) 



This is sufficient to duplicate the AS numerical calcula- 
tion. 

Figure 1 of Azbel and Soven& shows results from cal- 
culations made for three values of wavenumber k: 0.02, 
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FIG. 8: Localization length £ as a function of impurity con- 
centration c for two systems with different wavenumber k. 
Open symbols are calculations using the MK method, and 
solid lines are Eqn. 1271 The figure can be compared directly 
to Figure f of Azbel and Soven— 
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FIG. 9: The ratio of the measured localization £ to the cal- 
culated localization £ as a function of product of impurity 
concentration c and wavelength A. Filled symbols are calcu- 
lated solutions for systems having impurity cross section 0.2. 



1.5, 3.13. The MK method is used here to duplicate the 
numerical results for the two smaller values of fc, and the 
results are shown as open symbols in Fig. |S| The solid 
lines in Fig. [S] are the corresponding estimate of Eqn. 1271 
Two things should be noted explicitly: The AS definition 
of localization length corresponds to twice the localiza- 
tion length defined here. The definition of localization 
length in Azbel and Soven4 uses the Landauer— scaling 
parameter, while subsequent wor k 20,21 use the scaling of 
Anderson et ali£. The results shown in Fig. [S] use the 
latter scaling definitions. 



D. cA Effect 

As can be seen in Figs. ED and [7| the behavior of the 
total localization length diverges from the dilute limit £ • 
In Fig. EI the point at which £ begins to differ from £ is a 
function of the displacement wavelength. To more clearly 
demonstrate this effect, the ratio £/£ is shown in Fig. [5] 
as a function of the product cA. The data shown are those 
appearing in Fig. with the addition of those for A = 
512. Also shown in the figure is a dotted line denoting the 
long wavelength limiting behavior for A — > oo. For values 
of cA greater than 1, the dilute limit approximation does 
not hold, and the observed localization length is greater 
than the dilute limit estimation. Therefore, even though 
the effect is less dramatic for increasing A, there is a 
minimum effect, regardless of the wavelength. 



VI. ANHARMONIC CHAINS 

In an anharmonic system, phonon-phonon interactions 
will lead to the creation of displacements with varying 
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FIG. 10: The ratio of the measured localization £ to the di- 
lute limit localization length £ , normalized by the impurity 
concentration c, as a function of displacement wavelength A. 



wavelength. In time, very long wavelength displacements 
will be created. Because the impurity masses are fixed 
in time, the scattering cross section for the displacement 
waves will decrease as A~ 2 . In addition, because the im- 
purities are fixed in space, the product cA will increase 
over time. Based on the results from the last Section, one 
would expect that the ratio £/£ will diverge to infinity 
with increasing wavelength. 

Using Eqn. E| the ratio £/£ was calculated as a func- 
tion of displacement wavelength for systems having con- 
stant impurity mass and concentration. The results are 
shown in Fig. llUl for different values for m+. For all values 
of impurity concentration c, the ratio £ /£ asymptotes to 
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a constant at long wavelength: 

£ c 
lim = 1 + —— (m Q + m+) (29) 

Therefore, for anharmonic systems, the localization 
length is characterized, to within a constant, by the di- 
lute limit expression for £ given in Eqn. 1191 

The same relationship applies to both positive and neg- 
ative values for m+. For negative values of m + , in fact, 
in the limit m+ — * — 1, the behavior of the system is 
accurately characterized by £ D at long wavelengths. 

VII. CONCLUSION 

The localization length for harmonic chains having bi- 
nary disorder can be predicted accurately over a wide 
range of wavelengths, impurity cross section, and impu- 
rity concentration. The ingredient needed for this pre- 
diction is the cross section of a single scatterer, corrected 
for short wavelength displacements via the relative pulse 
propagation speed. The localization length over the en- 
tire impurity concentration range can then be estimated 
by approximating the system as a sum of two indepen- 
dent systems, each accounting for the behavior of the 
system at the two limits of impurity concentration. 

The general result applies to systems in which the im- 
purity mass is either larger or smaller that the original 
mass. Although the scattering cross section is symmet- 
ric about zero, with respect to the mass added to the 



background value, the localization length behavior dif- 
fers for negative and positive changes in mass having the 
same scattering cross section. This difference in behavior 
with respect to localization length, along with previous 
results showing differences in the rate of phonon-phonon 
interactions, are discussed in the context of numerical 
experiments on anharmonic systems. 

The general results also suggest that the localiza- 
tion length of long displacement wavelengths created 
by phonon-phonon interactions can be approximated, to 
within a constant, from dilute limit calculations results. 
For increasing concentration and constant impurity con- 
centration, the localization length, with respect to the 
dilute limit prediction, will eventually diverge toward in- 
finity. This deviation occurs for all wavelengths, and is a 
universal function of the impurity concentration and the 
displacement wavelength, for a constant cross section. By 
contrast, in a numerical experiment on an anharmonic 
system composed of fixed scatterers, the scattering cross 
section decreases with increasing wavelength. For these 
systems, the long wavelength behavior is, to within a con- 
stant, accurately predicted by dilute limit predictions. 
This conjecture will be confirmed in a future paper—. 
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